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We investigate the dynamical properties of Na 2 lrC> 3 . For five effective models proposed for 
Na 2 lr 03 , we numerically calculate dynamical structure factors (DSFs) with an exact diagonalization 
method. An effective model obtained from ab initio calculations explains inelastic neutron scatter¬ 
ing experiments adequately. We further calculate excitation modes based on linearized spin-wave 
theory. The spin-wave excitation of the effective models obtained by ab initio calculations disagrees 
with the low-lying excitation of DSFs. We attribute this discrepancy to the location of Na 2 lrC >3 in 
a parameter space close to the phase boundary with the Kitaev spin-liquid phase. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.70.Tj, 75.10.Kt 


I. INTRODUCTION 

Magnetic properties in 4 d and 5 d transition metal com¬ 
pounds have attracted much attention in condensed mat¬ 
ter physics. In some materials, such as Na2lr03, the 
energy scales of spin-orbit interactions, on-site coulomb 
interactions, and crystal fields compete with each other. 
The competition can produce unusual phases, including 
topological insulators (TJ and Kitaev spin liquids [ 2 j. 

In Na2lr03, Ir 4+ ions can be expressed as an isospin 
with a total angular momentum of 1/2 Q. We refer 
to this isospin as ’’spin” hereafter. IrC >6 octahedrons in 
Na2lrC>3 form a planar structure parallel to the ab plane, 
and Ir 4+ ions constitute a honeycomb lattice SI. In 
addition, IrOe octahedrons are connected by sharing the 
oxygen atoms on the edges, making the Ir-O-Ir bond an¬ 
gle nearly 90 °. This causes three kinds of anisotropic 
interactions between Ir 4+ ions depending on the bonding- 
path direction. Ir 4+ ions can also interact via direct over¬ 
lap of their orbitals. Thus, both Kitaev and Heisenberg 
interactions occur between Ir 4+ ions S, leading to the 
Kitaev Heisenberg model. 

Na2lr03 undergoes a magnetic phase transition to a 
zigzag antiferromagnetic order at Tn ~ 15 K a a. 
In a typical Kitaev-Heisenberg model, where ferromag¬ 
netic Kitaev and antiferromagnetic Heisenberg interac¬ 
tions are summed for nearest neighbor pairs, the zigzag 
order is not stabilized; thus, several models have been 
proposed to explain the zigzag ordering Hi- Some 
models S Hi have succeeded in explaining the tem¬ 
perature dependence of thermodynamic quantities, such 
as the specific heat and magnetic susceptibility. In dis¬ 
cussing interaction parameters, ab initio calculations are 
particularly powerful 1H3- However, there are consid¬ 
erable differences between the estimated parameters be¬ 
cause they are sensitive to the approximations used in 
the calculations. Therefore, there is still controversy sur¬ 
rounding suitable models for Na2lrC>3. 

The dynamical properties of Na2lr03 have been in¬ 
vestigated by inelastic neutron scattering (INS) exper¬ 
iments 0 , and a linearized spin-wave analysis has ex¬ 


plained the low-lying excitations observed in the experi¬ 
ments 0 , 7 ij. However, the proposed parameters differed 
between the studies. Choi et al. Q discussed the impor¬ 
tance of the long-range interactions, whereas Chaloupka 
et al. [ 7 ]] proposed the other scenario in which the signs of 
the Kitaev and Heisenberg terms play a key role. In ad¬ 
dition, the ab initio calculations indicated that Na2lrC>3 
is located close to the phase boundary with the Kitaev 
spin-liquid phase 00 - If the system is located close 
to the phase boundary, degenerate low-lying excitations 
from the magnetic frustration caused by the dominant 
Kitaev couplings may make conventional spin-wave the¬ 
ory invalid. Therefore, it is important to investigate the 
dynamical properties of Na2lrC>3 by using a method that 
does not depend on an approximation. 

In this paper, we numerically investigate the dynam¬ 
ical properties of Na2lrC>3 by an exact diagonalization 
method. We focus on dynamical structure factors (DSFs) 
of five effective models proposed for Na2lrC>3 
The DSFs provide magnetic excitations, which can be 
measured in INS experiments. We compare our numeri¬ 
cal results with the experimental results for the powder 
samples of this compound { 4 j and discuss the suitabil¬ 
ity of the models. To examine the low-lying excitations, 
we study excitation modes further by a linearized spin- 
wave analysis. The spin-wave excitations of the models 
obtained by the ab initio calculations disagree with the 
low-lying excitations of the DSFs. We attribute this dis¬ 
crepancy to the location of Na2lrC>3 close to the phase 
boundary with the Kitaev spin-liquid phase. Indeed, for 
an ab initio model [ 91 ], we confirm a double peak struc¬ 
ture in the specific heat, which can be a probe to observe 
the fractionalization of quantum spins in the Kitaev spin- 

liquid phase EH- 

The layout of this paper is as follows. In Sec. 
hd we introduce five effective models proposed for 
Na2lrC>3. We compute the DSFs by the numerical exact- 
diagonalization method to discuss the low-lying excita¬ 
tions of the five models. In Sec. m we discuss the low- 
lying excitations of the DSFs comparing with INS results 
for powder samples [ 1 ]. To capture the properties of the 
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FIG. 1: (Color online) (a) Honeycomb lattice model. Red, 
blue, and green lines denote Z- 1 X-, and y~ bonds, respec¬ 
tively. Dotted (dashed) lines represent the second (third) 
neighbor interactions. The colors of these lines are the same 
as those for the corresponding nearest neighbor lines. Dif¬ 
ferent lattice geometries for TV = 24 in (b) the 7 r/ 3 -rotational 
symmetric case (C 3 ), (c) 3 x 2 cells (3x2), and (d) 2x3 cells 
(2x3). The x, y, and 2 axes in (b) correspond to the orthog¬ 
onal axes for spin operators that are defined from the 5 d t 2 g - 
orbitals of Ir 4+ [g(] and the honeycomb plane is perpendicular 
to the (1,1,1) direction. Black bold dashed lines enclose the 
24 sites, and shaded rectangles indicate four sublattice unit 
cells. In (d), a schematic of the zigzag order is presented. For 
all geometries, we apply periodic boundary conditions and 
boundaries with the same symbol are connected. 


low-lying excitations, we calculate spin-wave excitations. 
We find that the spin-wave excitations fail to explain the 
low-lying excitations of the DSFs when the model is lo¬ 
cated nearby the Kitaev spin liquid phase. Finally, we 
conclude the discussion in Sec. m 


II. MODEL AND METHOD 
A. Effective models for Na 2 lr 03 


We consider a generalized form of Kitaev-Heisenberg 
models on a honeycomb lattice. The Hamiltonian is given 
as 


n = EE E ^Si-Sj 

p r p (ij)GTp 

+ EE E [Kr p SSSS + D rp (S l a S/ + h.c.)](l) 

p r p <ij)er p 


where T p is the bond indices depending on the direction 
( X , y , Z) of the p-th neighbor pairs and cr, / 3 , and 7 are 
indices for S = 1/2 SU( 2 )-spin components and take x, 
y , or z cyclicly. Each assignment of the indices depends 
on the bond direction; for example, 7 = z, a = x and 
/3 = y for the Z bond direction and so on. From the 
symmetry of the crystal structure for Na2lrC>3, the highly 
generalized form of the Hamiltonian reads 

^ = EE E E 

p r p (b)er p p.,v=x,y,z 


where the exchange coupling between i and j sites on the 
bond T p is given by a 3 x 3 matrix . In this paper, 
we focus on the five models EL 0-113 shown in Table Q] 
We can summarize three key factors in the five models: 
long-range interactions, trigonal distortions, and signs of 
the Kitaev interaction. In Models I and II, the inter¬ 
actions were evaluated by a spin-wave analysis so as to 
reproduce the low-lying excitations observed in the INS 
experiments and the temperature dependence of thermo¬ 
dynamic quantities 0 , 0 - In contrast, the interactions in 
Models III-V were estimated from the ab initio calcula¬ 
tions [IHI3. In Models III and V, interactions were es¬ 
timated from the density-functional-theory calculations. 
Note that there are technical differences between Models 
III and V in evaluating the tight-binding models and the 
consequent effective spin model. In contrast, in Ref. Q, 
Katukuri and co-workers employed ab initio techniques 
from wave-function-based quantum chemistry and pro¬ 
posed several parameter sets for the coupling constants. 
We adopt the nearest-neighbor interactions that were 
used in the right phase diagram in Fig. 2 of Ref. Q as 
Model IV. For the second and third neighbor Heisenberg 
interactions, we adopt the middle values between the pro¬ 
posed range for explaining the experimental Curie-Weiss 
temperature [§]. Details of the interaction parameters for 
the five models are summarized in Tables eh and EU 


B. Dynamical structure factors 


For the five models presented in the previous subsec¬ 
tion, we calculate the DSF, S(Qi,ix), for the system size 
TV = 24 with three different lattice geometries shown in 
Figs. Gib) Old). The DSF at a zero temperature is de¬ 
fined as 


S^(Q, U ) = -~ lim lm( 0 o| 

tt £->-+0 to + Eg + le — ft 


I^o), ( 3 ) 


where </>o is the ground state of H with the energy Eq and 
Sq = exp(— iQ ■ r). </o and Eg are calculated 

by the Lanczos method, and then S IJV ( Q. ui) is obtained 
by a continued fraction expansion fl 3 l . fTIj . 

In general, when the off-diagonal elements of are 
non-zero, the DFS S flv (Q,uj) is allowed to have non¬ 
zero off-diagonal elements. The contribution from the 
off-diagonal elements in the DFSs is expected to be pro¬ 
portional to the off-diagonal spin correlation. The am¬ 
plitude of the spin correlation is also proportional to the 
absolute values of the matrix element of J^ v ■ In the 

1 P 

five models presented in the previous section, the sum of 
all diagonal elements in is larger than the remain¬ 
ing each elements. Therefore, we consider that the sum 
of diagonal elements, namely S(Q,lu) = ^ ^^(Q, w), 
mainly contributes to the scattering intensity. 

The magnetic excitations of Na2lrC>3 have been in¬ 
vestigated by INS experiments for the powder sam¬ 
ples 0 . Therefore, the scattering intensity observed 
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TABLE I: Interactions and effects of trigonal distortions in the five proposed models. l y/’ —indicates that the corresponding 
term is included (not included). V indicates that it was present, but the value was not estimated. F (AF) denotes ferromagnetic 
(antiferromagnetic) Kitaev interactions. NN (NNN) denotes next-(next-) nearest-neighbor interactions. LSW: linearized spin- 
wave theory. DFT: density-functional theory. QC: wave-function based quantum chemistry. 


Model 

Long-range interaction Trigonal distortion Kitaev interaction 

Method 

Model I: [4| 


— - 

LSW 

Model II: [7] 


AF 

LSW 

Model III: [9] 

V 

V F 

ab initio [DFT] 

Model IV: [8] 

* 

V F 

ab initio [QC] 

Model V: [10] 

V 

F for NN and AF for NNN 

ab initio [DFT] 


TABLE II: Coupling constants except for Model III. The pos¬ 
itive (negative) number corresponds to the antiferromagnetic 
(ferromagnetic) interactions. For Models I, II, IV, and V, 
the Hamiltonian can be expressed by the expression © with 
D 2nd = D 3ld = 0. The energy unit is meV. 


^ylst ^-lst 

D ist 

j2nd 


J' ild 

K'drd 

Model I: [4] 4.17 0 

0 

3.25 

0 

3.75 

0 

Model II: [7] -4.0 21.0 

0 

0 

0 

0 

0 

Model IV: [24] 3 -17.5 

-1 

4.5 

0 

4.5 

0 

Model V: [10] 5.8 -14.8 

0 

-4.4 

7.9 

0 

0 


in the experiments is averaged with respect to the 
wave vectors and the scattering directions. In or¬ 
der to compare the experimental results with numeri¬ 
cal results, we consider the averaged intensity defined 
as l(\Q\,w) = Yj\ Q\=Qi S (Qi,u), where S{Qi,ui) = 
y z w) and the sum for Q, runs all pos¬ 

sible wave vectors inside the first Brillouin zone. Thus 
the averaged intensity, I{\Q\,w), is characterized as a 
function of the distance from T point. 


III. RESULTS 

A. Ground states of five models 

Figure [2] shows the results for the static structure 
factors (SSFs) for the five models. The longitudinal 
and transverse elements of the SSFs are defined as 
S Z (Q) = f 0 °° S zz {Q,uj)dLo and S‘(Q) - J™(S xx {Q, u) + 
S vv (Q,uj))du}, respectively. Although the amplitude of 
the SSFs depends on the lattice geometry, the largest 
peak appears at the Y and M points in Models I -V. The 
results indicate that the stable ground state of the five 
models is a zigzag order. This is consistent with the ex¬ 
perimental results SI- 

Because we expect that the results of the tt/3- 
rotational-symmetric lattice, which is labeled C 3 in Fig. 
[ 2 ] well describe the properties at the thermodynamic 
limit, we focus on the results for the C 3 case. The ground 
states, except for in Model III, are characterized by a 
zigzag spin configuration, although the ’’type” discussed 


below cannot be assigned because of the C 3 three-fold 
rotational symmetry. In Model III, the largest peak ap¬ 
pears at the Y point in the transverse component. This 
means that the spin correlation on the Z bond is anti¬ 
ferromagnetic, and thus the ground state exhibits Z-type 
zigzag order [9j. Moreover each magnetic moment of the 
ground state points to a direction in the plane perpen¬ 
dicular to the z axis [Fig. ©b)]. 

B. Powder averaged results 

Figure [3] shows the numerical results for the aver¬ 
aged intensity, I{\Q\,w). We compare X(|Q|,w) of the 
five models with the INS results for the powder samples 
shown in Fig. 3 of Ref. Q. The INS experiments Q show 
that there are three characteristic features for T < T^: 
(i) a sharp lower boundary for the scattering intensity 
below 4 meV between the T point and Y point (except 
for IQI ~ 0 and ui < 2 meV, where the reported data 
have been lacked), (ii) strong scattering peaks at close to 
4 meV between the Y point and T* point, and (iii) no 
strong intensities at the T* point Q~5f. 

First, we discuss the excitation boundary. The excita¬ 
tion boundary begins at 4 meV nearly midway between 
the r point and Y point, and the boundary energy de¬ 
creases towards the Y point [dj. The low-lying excita¬ 
tion energy in the related region for Models I-IV draws 
the convex curvature and it agrees with the experimen¬ 
tal results. In contrast, the low-lying excitation energy 
in Model V also increases from the Y point. However, 
it exhibits a strong lattice-geometry dependence and the 
convex curvature can not be observed. The maximum en¬ 
ergies of such curvature for Models I -IV locate at nearly 
midway between the T point and Y point and are about 
3.8 meV, 6.4 meV, 2.6meV, and 5.3 meV, respectively. 
Therefore, the agreement is excellent for Model I. Since 
Model I was estimated to explain the excitation bound¬ 
ary by the linearized spin-wave theory, the agreement is 
quite natural. Despite the good agreement, Model I is 
inadequate for describing Na 2 lrC> 3 , because it does not 
include the Kitaev interactions. It is unlikely that there 
is no Kitaev term, when we consider the interaction path 
between Ir 4+ ions in Na 2 lr 03 . 

Secondly, we focus on the peaks below T = 7 n de- 
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TABLE III: Coupling constants for Model III. Model III is expressed by the Hamiltonian ©• The energy unit is meV. 
means that the coupling is not included. 


X 

T 1st 

J X 

y 

2 

rlst 

J y 
x y 

z 

T 1st 

j z 
x y 

2 

j'Z nd 

J X 

x y z 

T2nd 

J y 

x y z 

X 

j2nd 

Jz 

y z 

X 

T3rd 

J X 

V 

2 

X 

7-3 rd 

J y 

y 

2 

X 

T3rd 

J Z 

V 

2 

x -23.9 

-3.1 

-8.4 

2.0 -3.1 

1.8 

4.4 -0.4 

1.1 



-0.8 

1.0 -1.4 

1.7 

0 

0 

1.7 

0 

0 

1.7 

0 

0 

y -3.1 

3.2 

1.8 

-3.1 -23.9 

-8.4 

-0.4 4.4 

1.1 

- 

- 

1.0 

-0.8 -1.4 

0 

1.7 

0 

0 

1.7 

0 

0 

1.7 

0 

s -8.4 

1.8 

2.0 

1.8 -8.4 

3.2 

1.1 1.1 

-30.7 



-1.4 

-1.4 -1.2 

0 

0 

1.7 

0 

0 

1.7 

0 

0 

1.7 



r x r x r xr xr x 


a 


a 


a 


a 


a 


FIG. 2: (Color online) (a) Four sublattice unit cells and schematic spin configuration of the Z-type collinear zigzag order. Solid 
(open) circles represent up (down) spins and the dotted rectangle corresponds to the shaded ones in Fig. [T] (b) Reciprocal 
vectors a and /3 for the four sublattice unit cells in (a), (c)-(l) SSFs for Models I-V for N = 24 with three lattice geometries in 
Figs. ©ty-IHd). The area of each circle is proportional to the amplitude. Upper (lower) panels are results for the longitudinal 
(transverse) component. 


scribed in the experimental feature (ii). When we com¬ 
pare the experimental results for T < Tjy with those 
for T > Tn, we observe scattering intensities derived 
from the magnetic ordering around 4 meV between the 
Y point and T* point. Our numerical results reproduce 
these scattering intensities, except for Models II and IV. 
However, in Models I and V, relatively large peaks ap¬ 
pear around ~ 3.8 meV at the T* point. The presence 
of such peak contradicts the experimental feature (iii). 
Thus, Model III is also considered as a good candidate 
for the model that explains the experimental features (i), 
(ii), and (iii). Therefore, among the five models, Model 


III is the most suitable for explaining the INS experi¬ 
ments for Na 2 lr 03 . However, the low-lying excitation of 
the DSFs in Model III appears slightly lower than the 
experimental results [dj. Thus, further examination of 
the second and third neighbor interactions is desirable to 
improve the accuracy of the theoretical prediction. 


C. Low-lying excitations 

Next, we discuss characteristics of the low-lying excita¬ 
tions of the DSFs, S(Qi,u>). We compare the DSFs with 
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FIG. 3: (color online) Averaged intensity, T(\Q\,uj), for N = 24 with three lattice geometries. The area of each circle is 
proportional to the scattering intensity. All results are normalized by the largest value for each data set. Dotted lines are the 
low-energy cutoff at about 2 meV in INS experiments Q. The horizontal axis corresponds to distance from the T point. The 
labels on the horizontal axis denote the corresponding distance from the T point to the label [li¬ 


the spin-wave excitation modes for the four sublattice 
unit cells [Figs. [H[a)-[IIe)]. In the linearized spin-wave 
calculations for Models I, II, and V, we assume the Z- 
type collinear zigzag configuration shown in Fig. [2](a). 
We examine the initial state of the spin-wave analysis for 
Model V, namely the ground state for the classical spin 
model. Although Model V includes the long-range in¬ 
teractions expressed by the Kitaev- Heisenberg term, the 
ground state analysis for classical spins indicates that 
the collinear zigzag order is still favored. For Models 
III and IV, we obtain a collinear zigzag spin configura¬ 
tion for the ground state from the classical spins analysis. 
However the zigzag order parameter, which is defined as 
M = Si + S 2 — S 3 — S 4 for four spins in the unit cell in 
Fig. da), indicates that they are in a different direction 
from those for Models I, II, and V. This tilting of the 
order parameter results in the spin gap at the M point. 

We turn our attention to the results for Model III. The 
spin-wave excitations of Model III appear in the higher 
energy region above lu/Jnn ~ 1 with an energy gap at 
the M point. The low-lying excitation energies of the 
spin-wave modes at other Q s are a few-times higher than 
those of the DSFs. Thus, in Model III, the conventional 
spin-wave theory clearly fails to describe the low-lying 
excitations of Na 2 lrC> 3 . 

A similar breakdown of the conventional linearized 
spin-wave picture for the low-lying excitations of the 
DSFs is also observed in Models IV and V, where the pa¬ 
rameters are also estimated from ab initio calculations. 


In Model V, except for the M point, the low-lying exci¬ 
tation of the spin-wave mode at each Q also appears in 
the high-energy region. We here comment on the high- 
energy excitations in the DSFs. We observe some poles 
in the DSFs at similar energies to the spin-wave excita¬ 
tions, although those are not shown in Fig. [4] owing to 
the large discrepancy in the energy scale. The observed 
breakdown of the spin-wave excitations indicates that the 
low-lying excitations of the DSFs in Models III-V are ap¬ 
parently different from the free magnon excitations and 
can be attributed to precursors to deconfined Majorana 
excitations, which are composed of itinerant Majorana 
fermions and Z 2 gauge fields [H, [23[. 

IV. DISCUSSIONS 

The spin-wave picture fails to explain the low-lying ex¬ 
citations in ab initio models for Na 2 lrC> 3 . This becomes 
significant when the model is close to the boundary of the 
spin-liquid phase, as for Model III [g]. Recently, a similar 
breakdown of the spin-wave picture has been discussed 
in INS measurements for a potential Kitaev material, a- 
RUCI 3 [l 6 |. This material is near the Kitaev spin-liquid 
phase boundary [TtI l22j| . Close to the Kitaev spin-liquid 
phase boundary, we expect that a precursor of the frac- 
tionalization of quantum spins into itinerant Majorana 
fermions and Z 2 gauge fields [lH, |23j] will be observed 
and may break the spin-wave picture for low-lying exci- 
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FIG. 4: (Color online) DSFs, S(Qi,u), for N = 24. Results for different lattice geometries are included. The area of each circle 
corresponds to the logarithmic amplitude of the scattering intensity. The solid curves are dispersion curves estimated from the 
linearized spin-wave analysis. The excitation energy is normalized by the largest Heisenberg interaction | Jnn\ in each Model; 
| J iviv j=4.17 meV, 4.0 meV, 4.4 meV, 3.0 meV, and 5.8 meV for Models I-V, respectively. The horizontal axis in (a)-(e) runs 
along the arrows shown in Fig. [2] (b). The excitation energy is normalized, (f) Temperature dependence of specific heat for 
Model III numerically obtained by a diagonalization method. Red (blue) symbols are the results for IV = 16 (8). 


tations. The double peak structure in the specific heat 
can be used as a probe to observe the emergence of the 
fractionalization in the Kitaev spin liquid for the current 
two-dimensional models (l3| . In Model III, we confirm 
the double peak structure in the temperature dependence 
of the specific heat [Fig. QJf)]- Thus, the present results 
demonstrate that Na 2 lr 03 is located close to the Kitaev 
spin-liquid phase boundary, if Model III is the most suit¬ 
able. 

The trigonal distortion is also important for the ‘dis¬ 
tance’ from the Kitaev spin-liquid phase in Model III Q. 
We believe that systematic comparisons with the energy 
scale of the low-lying excitation of the DSFs, the peak 
separation of the specific heat, and the trigonal interac¬ 
tion may clarify discussions of the Majorana physics in 
experimental observations. 
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